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We study the relaxation of a single electron spin in a circular gate-tunbable quantum dot in gapped 
graphene. Direct coupling of the electron spin to out-of-plane phonons via the intrinsic spin-orbit 
coupling leads to a relaxation time Ti which is independent of the B-field at low fields. We also 
find that Rashba spin-orbit induced admixture of opposite spin states in combination with the 
emission of in-plane phonons provides various further relaxation channels via deformation potential 
and bond-length change. In the absence of valley mixing, spin relaxation takes place within each 
valley separately and thus time-reversal symmetry is effectively broken, thus inhibiting the van 
Vleck cancellation at B = known from GaAs quantum dots. Both the absence of the van Vleck 
cancellation as well as the out-of-plane phonons lead to a behavior of the spin relaxation rate at low 
magnetic fields which is markedly different from the known results for GaAs. For low B-fields, we 
find that the rate is constant in B and then crosses over to oc or oc at higher fields. 



Introduction — The electronic spin degree of freedom is 
under intense investigation as a possible implementation 
of a qubit [T]. While the feasibility of all required oper- 
ations has been experimentally demonstrated for GaAs 
quantum dots (QDs) [2, the decoherence caused by the 
surrounding nuclear spins in the host material remains 
challenging. Regarding the use of the electron spin as 
a qubit in quantum computation devices, spin decoher- 
ence and relaxation are limiting factors. In general, a 
necessary condition for a working qubit is that the time 
required to perform an operation is significantly shorter 
than the decoherence and relaxation times. Motivated by 
this, the implementation of qubits in QDs in graphene 
has been proposed [S]- Graphene consisting of natural 
carbon comprises 99% of the carbon isotope ^^C which 
has no nuclear spin, hence the hyperfine interaction is ex- 
pected to play only a minor role. Furthermore, spin-orbit 
interaction (SOI) in graphene is expected to be relatively 
weak and therefore long decoherence times are expected. 
However, for spins localized in QDs in carbon nanotubes, 
SOI has turned out to be unexpectedly strong [H [5] due 
to curvature-induced effects. It has also been shown the- 
oretically that van Hove singluarities in the phonon den- 
sity of states in one dimension can lead to strong varia- 
tions in the spin relaxation rate [6] . It is therefore impor- 
tant to investigate the spin relaxation time in graphene 
QDs. The form of the SOI in graphene, both intrinsic and 
Rashba type, is known [7], and there are several works 
on its strength depending on various parameters such as 
curvature or electric field [H E]. There have also been 
experimental [TOl and theoretical pU] studies on spin re- 
laxation of extended states in graphene. 

In this paper we determine theoretically the spin relax- 
ation time Ti for an electron confined to a circular QD in 
gapped graphene as a function of the external magnetic 
field B. It has been predicted previously that such QDs 
can be formed with electrostatic confinement in either 
single-layer graphene with a substrate-induced band gap 




FIG. 1: (Color online) a) The two states of a spin qubit 
(blue solid arrows) reside in the same valley, as opposed to a 
Kramers qubit (empty red arrows) , formed by a Kramers pair 
related by time-reversal symmetry (T). While in single- valley 
semiconductors such as GaAs these two cases are identical, in 
graphene the Kramers qubit involves states in different valleys 
{K and K'). b) The B-field orientation is given with spherical 
coordinates 9 and (pB relative to the normal to the graphene 
plane. The propagation direction of the emitted phonon (red 
wavy arrow) is described by the angle (j>q. 



or bilayer graphene with an electrically controlled gap 
jl2| . At B = 0, the states in these QDs have a two- 
fold valley degeneracy which can be lifted in a perpen- 
dicular magnetic field. Being a centro-symmetric crystal, 
phonons in graphene do not couple piezo-electrically, thus 
leaving three possible electron-phonon coupling (EPC) 
mechanisms: deformation potential, bond length change, 
and direct spin-phonon coupling. From these EPC mech- 
anisms, we derive two spin relaxation mechanisms. One 
such mechanism involves the admixture of states of oppo- 
site spin and excited orbitals into the dot eigenstates due 
to SOI, in combination with energy relaxation via phonon 
emission [131 [T3] . It turns out that to lowest order in the 
EPC, this only involves in-plane phonons coupled via the 
deformation potential and bond-length change. The sec- 
ond mechanism directly couples the spin to out-of-plane 
phonons via curvature induced SOI. For comparison, in a 
parabohc GaAs QD, a strong dependence oc has been 
predicted for both mechanisms fl3l. Relaxation times in 
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the millisecond range at a field B = IT have been pre- 
dicted and even longer Ti exceeding one second have been 
experimentally verified |15j . The prediction for graphene 
QDs looks markedly difi'erent because of the absence of 
the van Vleck cancellation for spin qubits in a single val- 
ley as opposed to "Kramers qubits" (see Fig. [I^) , as well 
as the absence of piezo phonons and the two-dimensional 
nature of phonons. 

Model — To study spin relaxation in a circular and gate- 
tunable QD in single-layer graphene, we assume the host 
graphene layer to be sufficiently large to ensure that the 
edges do not induce inter-valley mixing. The QD can be 
described with the Hamiltonian pj], 



Ho = vf{p + eA_L) • cr + ^g^lB'B ■s + U{r) + rAa,, (1) 



where the first term is the well-known Dirac Hamiltonian 
for graphene in the presence of a vector potential 
Aj_ with B_L = V X A_L = (0, 0, B cos 6) being the per- 
pendicular component of an arbitrarily oriented B-field 
(Fig. [ij)). The second and third terms describe the Zee- 
man coupling of the electron spin to the total B-field and 
the smooth and circularly symmetric confinement poten- 
tial U{r). The last term opens a band gap 2 A which can 
e.g. arise due to the influence of the substrate [T?! [T8]. 
Here, r = ±1 denote the K and K' valleys. In the ab- 
sence of valley-scattering, we can restrict ourselves to a 
single valley, e.g. r = -1-1. Weak inter-valley coupling 
can arise from atomic defects or boundaries |16| . or from 
the hyperfine interaction with the remaining ^^C atoms 

m- 

The eigenstates |n, s)'^"^ of Hq in Eq. ^ with energy 
En + sgfiBB/2 are simultaneously eigenstates of the to- 
tal angular momentum j G 'Z + ^, i.e. orbital quantum 
number and pseudo-spin, with spinor wavefunctions (in 
the K valley). 



Xa (r) 



(2) 



The spinor components Xa^i^) can be given in closed 
form for a step-like potential U{r) — Uod{r—R) \12\, how- 
ever the eigenenergies En have to be evaluated numeri- 
cally. Each eigenstate is characterized by a pair (n, s) 
where s =t,J,= ±1 is the spin and where n = (j^, j) has 
a radial and angular momentum part. 

In-Plane Phonons — In order to study processes based 
on the admixture mechanism, we begin with the Hamil- 
tonian H = Hf) + Hso + ^^EPC , where Hq describes the 
graphene QD without SOI as explained above, ffso de- 
scribes the SOI, and Hepc describes EPC. The effect of 
the SOI is to weakly mix the eigenstates Eq. ([2]). In this 
manner, e.g. the QD ground state, say |n = (0, 1/2), t)'"-* 
acquires components of the excited states with 



n' = (v'^j') 7^ n and opposite spin, to first order in -ffsOj 

n'^n ^ 

and similarly for |n i). With this admixed state the spin- 
conserving EPC can cause spin relaxation. 
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■|ffEPc|n;)^(i^EPc)It = 

(^^So)It' (^^EPC)„'„ , (^EPC),: 



{Hso)]^n 



En — En' — ^gfisB En — En' + ^gfisB 



(4) 



For sufficiently small B-fields this can be expanded 
around B — 0. In the case of GaAs, the expres- 
sion Eq. Q vanishes for B = due to the symmetry 

{Hso)ntn = ~ i^So)mn ^^'^ {HEPc)nm = {HEPc)mn- 

This van Vleck cancelation |T31 [T3] is one of the rea- 
sons for the high power of B that appears in the spin 
relaxation rate in GaAs QDs and can be traced back to 
the time-reversal invariance of H and its eigenstates, i.e., 
the fact that both SOI and EPC preserve time-reversal 
invariance. In particular, the spin relaxation takes place 
from one state, say \n'[)^ to its partner jnj,) within a 
Kramers pair, which are linked by time reversal. 

In our case, the states jnf) and \nl) lie in the same 
valley and therefore do not form a Kramers pair (see 
Fig. [T^). The time-reversed partner of jnt) is \nl)', 
where the prime denotes the opposite valley. Since nei- 
ther the EPC nor the SOI lead to inter-valley mixing, 
spin-relaxation is effectively constrained to a single val- 
ley. Therefore the selection of spin qubit states within 
the same valley breaks time-reversal symmetry and leads 
to the absence of the van Vleck cancelation. We now 
proceed to the evaluation of the matrix elements of the 
SOI and the EPC in Eq. Q in order to calculate the spin 
relaxation rate. 

We divide the SOI Hamiltonian into its intrinsic and 
Rahsba terms [7], 

Hso ^ Hi + Hji^ AjTCT^s^ -I- A^,{TaxSy - ayS^), (5) 

where ai and Si denote the Pauli matrices acting on 
the pseudo-spin and real spin. We use a spin quanti- 
zation axis aligned with the external B-field (see Fig.[T]D) 
and corresponding spinors | fs) and | Ib) and obtain 
= (tBk.lis) = cos2 I - sin^ f and fy = (fs 



sin^l; 



First we consider 



\Sy\iB) ^ -i(cOS^ 2 

Hn and calculate its matrix elements with states \n f^) 
and \n' Ib)- The two spin states we use are orthogonal, 
i.e. (ts I Ib) = but they are not Sz eigenstates. In 
principle this allows both the intrinsic and Rashba SOI 
to provide a relaxation channel in the admixture mech- 
anism. However, due to the circular symmetry of the 
dot, selection rules for j apply. In the case of Hji only 
dipole transitions (|j — j'\ = 1) are allowed, whereas the 
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matrix element of Hi gives rise to selection rules j = j' 
which turns out to be incompatible with the selection 
rule \j - j'l = 1 for the EPC. 

The matrix element of Hfi can be written as {Hji)]^^, ~ 
2^A^[/n<5,y+iiV„^„^ + 5,,'-iiV„^f ) - in5,r+iN^S - 
5,,v_iiV^^)], where N^^ = /drrxlxS'- The matrix 
element {Hjif^n' ^® neither symmetric nor antisymmetric 
in contrast to the case of GaAs where an antisymmetry 
leads to van Vleck cancelation. 

We consider two different EPC mechanisms which cor- 
respond to different changes in the lattice induced by 
phonons. The deformation potential is caused by an area 
change of the unit cell, whereas the bond-length change 
mechanism corresponds to a modified hopping propabil- 
ity [201 121j . Because we work in the low-energy regime, 
we only consider acoustic phonons. In principle there are 
six possible relaxation channels: (i) longitudinal acoustic 
(LA) , transversal acoustic (TA) , transversal out-of-plane 
(ZA) phonons, and (ii) deformation potential [gi) and 
bond-length change (172) mechanisms. In lowest order in 
the atomic displacement, the EPC has the form ^01 [21] 



3202 SlOi 



with ai = i and 02 



(6) 

'1 for LA phonons, and 02 = 
g2i0q and ai — Q for TA phonons, and A the area of 
the graphene sheet. The vanishing of ai is due to the 
fact that in the regime of linear atomic displacements the 
coupling of the TA mode is a two-phonon process. Here, 
we restrict our considerations to one-phonon processes. 
For a B-field of i? = IT and a sound velocity of s = 
2 X 10'* m/s [22], we obtain from g^isB = hsq a phonon 
wavelength of A w 300 nm which is an order of magnitude 
larger than a typical QD size of 25 nm [7, thus justifying 
the use of the dipole approximation for typical laboratory 
fields. 

For the matrix element for LA phonon coupling 
via the deformation potential we find (-ffEPc),™' ~ 



31-^ 



'+ie' 



(Sjj'.ie**") with 



Jdr r^xTx'X+X'h*X% 



The dependence 
g disappears upon 



on the phonon emission angle 
summation over final states. For the TA phonons we 
find that the coupling via the deformation potential is a 
two-phonon process which will not be discussed here. 

The bond-length change mechanism leads to similar 
results for both LA and TA phonons, (-ffEPc)„„' — 

-Dla = -i27rg2/V^psLA and L'ta = '^''^92/ \/ Aps^A, and 
where the plus (minus) sign corresponds to LA (TA). 
In linear order in the atomic displacement the ZA mode 
is decoupled from the other modes. The Hamiltonian 
Eq. ^ cannot account for a coupling to the out-of-plane 
mode. 

With the matrix elements derived above, we can write 





V = 1 


1^ = 2 


u = 3 


j = -0.5 

3 = L5 


1.1 X lO'* 
1.6 X 10" 


2.6 X 10"^ 
1.2 X 10^ 


1.3 X 10"^ 
9.3 X 10"^ 



TABLE I: Individual relaxation rates in units of s^^ for LA 
phonons via the deformation potential at B = 1 T. For higher 
quantum numbers v, the rate decreases quickly. 



the transition rates using Fermi's golden rule as l/Ti 
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2nAj 



(27r)2 



)n 
nn 



^(■s? — gfJ-sB). For all 



(-ffEPc)f 

mechanisms we find the same dependence on the ori- 
entation of the B-field, f{9) = cos4(6l/2) + sin*(6l/2) = 
(3 + cos(26'))/4. We find for the relaxation rate from the 
deformation potential 



ai 



(7) 



while for the bond-length change mechanism, we have 



92 



^LA.TA 



n' 



(%-+i(A^„^,^?)^ + %,_i(A^„^f)2) 



, (8) 



with Rnn' — {En — En')^^- Foi numerical evaluation, 
we assume a QD size of i? = 25 nm and A = 106 where 
6 = v/R is the average level distance. The depth of the 
quantum well is also set to Uq = 106. The Rashba SOI 
constant can be adjusted by an external electric field [9] 
or by using different types of substrates. We chose a 
value of A^ = 48 /ieV to calculate the relaxation times 
displayed in Fig. [2] For the EPC constants we assume 
gi — 30 eV and 32 = 1.5 eV [ID]. We use as sound veloci- 
ties SLA = 1.95 X 10'' m/s and sta = 1-22 x lO^m/s [32]. 
The overlap integrals and M„„/ are calculated nu- 

merically. The sum over n' runs over all states, including 
the continuum. As shown in Table |T] the contributions 
from higher levels vanish quickly so that we only take 
the first three levels into account. The relaxation rate 
Ti = l/F is plotted in Fig. [2) 

Direct Spin-Phonon Coupling — In flat graphene the 
acoustic phonons with perpendicular (ZA) polarization 
are decoupled from the in-plane modes (LA,TA). We 
extend the SOI Hamiltonian Eq. ([5| for the case of 
a graphene layer which is curved due to ZA phonons. 
For displacements much smaller than the wavelength the 
normal vector of the graphene plane can be written as 
n(z) « z -I- Vuz(a;,y) where Uz{x,y) is the displacement 
held representing the ZA-phonons. Rotating the spin ma- 
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FIG. 2: (color online) Log-log plot of the spin relaxation time 
Ti as a function of an external B-field perpendicular to the 
plane {6 = 0) defined by the graphene sheet. The radius of 
the dot is _R = 25 nm, both energy gap and depth of the dot 
are 260 meV. The individual relaxation channels are the cou- 
pling to LA in-plane phonons via deformation potential {gi, 
red dotted line) and the coupling to LA and TA phonons via 
bond-length change {gi, blue dashed and green dot-dashed 
lines), as well as the direct coupling to the out-of-plane (ZA) 
phonons (purple, long-dashed line). The black solid line rep- 
resents the sum of all four processes. Inset: Dependence of 
the relaxation rate on the inclination angle 9 of the B-field. 



trices into the local frame determined by the normal vec- 
tor n(z) we obtain in linear order in u{z) a generalized 
SOI Hamiltonian Hso = Hi + Hjx with 



Hi = H 



(0) 



Hk - H^ 



(0) 



(9) 
(10) 



where h['^^ and h'^^ are the SOI Hamiltonians for 
flat graphene given in Eq. ([5]). We evaluate these 
expressions for transverse out-of-plane (ZA) phonons, 
with a quadratic dispersion relation = ^q^ where 
/i = \/~kJp with K = 1.1 eV the bending rigidity and 
p = 7.5 ■ lO^'^'Icg/m^ the mass area density [22l |24] . 
The EPC Hamiltonian is then obtained by substituting 
the displacement operator for the ZAphonons — 
yi/ApWq (e^'i-''6t + e-*'i-'"6) into Eqs. ^ and ([lo|. For 
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the intrinsic SOI we obtain the matrix element [Hi) 
iAi./l/A^{qM \s^\ i) +%(t \sy\ i)) (n|a,e'q--|n). 
When evaluating the orbital matrix element only the 
lowest order in the dipole approximation contributes. 
All higher orders contain a factor c>c e^'^i which averages 
to zero when the integration over 0q is carried out. 

Finally, Fermi's Golden Rule is used to find the relax- 
ation rate 



^ZA 



dr r 



\Xb\ 



(11) 



which is independent of B. The Matrix element itself 
depend only weakly on B. For the numerical evaluation 



we use Ai = 12/zeV [13] and sza = 1-59 x lO^m/s 
The same calculation for the Rashba SOI yields vanishing 
matrix elements and therefore no additional contribution. 
In some cases, boundary conditions may lead to a linear 
dispersion relation for the ZA-phonons. We find that in 
this case the contribution due to ZA-phonons is negligible 
compared to the in-plane phonon contributions. 

Conclusion — We have calculated the electron spin re- 
laxation time Ti in a gate-tunable graphene QD aris- 
ing from the combination of SOI and EPC. We have re- 
stricted ourselves to the zero-temperature case, i.e. pure 
phonon emission which is realistic at 0.1 T and 100 mK 
and higher temperatures for larger fields. We have taken 
into account two mechanisms: Admixture mechanism 
and direct spin-phonon coupling. Due to selection rules 
in a circular QD, the admixture mechanism only leads 
to spin relaxation in combination with the Rashba SOI. 
The deformation potential EPC with LA phonons leads 
to a spin relaxation rate scaling as B'^ (Fig. [2| , while the 
bond length change EPC with both LA and TA phonons 
results in B^ dependencies. The relatively low powers 
compared to GaAs QDs can be traced back to the ab- 
sence of the van Vleck cancelation, in combination with 
the 2D phonon density of states. The direct coupling of 
electronic spins to ZA phonons only leads to spin relax- 
ation in combination with the intrinsic SOI whose rate 
does not depend on the applied B-field (in lowest order) 
and thus leads to a B-field dependence at low fields which 
is markedly different from that in GaAs QD. 
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